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i—i Abstract: 
> 

BSLP is a two-dimensional dynamic model of interacting portfolio-level loss and loss in- 
tensity processes. It is constructed as a Markovian, short-rate intensity model, which 
facilitates fast lattice methods for pricing various portfolio credit derivatives such as 
tranche options, forward-starting tranches, leveraged super-senior tranches etc. A semi- 
parametric model specification is used to achieve near perfect calibration to any set of 
consistent portfolio tranche quotes. The one-dimensional local intensity model obtained 
in the zero volatility limit of the stochastic intensity is useful in its own right for pricing 
^ non-standard index tranches by arbitrage-free interpolation. 
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1 Introduction 



A large class of portfolio credit derivatives can be viewed as derivatives referencing the 
cumulative portfolio losses. Furthermore, we can distinguish between two classes of such 
derivatives. For single-period (vanilla) instruments such as synthetic CDO tranches all 
that is needed for pricing is the set of marginal portfolio loss distributions at different 
time horizons. More exotic, multi-period instruments such as tranche options require 
knowledge of the future distributions of the mark-to-market (MTM) of the tranche as 
well as the portfolio losses. This is equivalent to having a model for the evolution of the 
portfolio losses, from which the MTM of the tranche can then be derived. 

If we hedge such multi-period products only using tranches on the underlying portfolio, 
then the challenge of modeling single-name dynamics does not arise at all. Therefore, 
under these conditions both the pricing and risk management of portfolio credit derivatives 
can be done within a model that specifies the portfolio-level dynamics of cumulative losses 
but leaves the single name dynamics unspecified (or specified at a later stage). 

The fact that powerful and flexible pricing models can be constructed within such 
a "top-down" framework^ was first recognized in work by Giesecke and Goldberg pQ, 
Schonbucher [2], Sidenius, Piterbarg and Andersen (SPA) [3], and Bennani [I]. More 
recent publications along this strand include Brigo et al. [5], Chapovsky et al. [6], Errais 
et al. [7j, Cont et al. [8], Ding et al. [9], and de Kock et al. [TP] . Furthermore, while this 
paper was under preparation we learned of a recent paper by Lopatin and Misirpashaev 
[TT] who have independently suggested a model similar to the one developed in this paper. 
We will comment on similarities and differences between their approach and ours in due 
course. 

In this paper, we present the Bivariate Spread-Loss Portfolio (BSLP) model - a dy- 
namic Markovian model for correlated portfolio loss and loss intensity processes. In our 
model the portfolio loss process follows a Markov chain whose generator is driven by a 
stochastic intensity (so that the generator itself becomes stochastic). The intensity is 
given by a diffusion process which can incorporate default-induced jumps. The fact that 
the driving intensity and the loss process are mutually dependent means that our frame- 
work is more general than the more standard doubly stochastic one which only allows for 
a one-sided dependence. 

The portfolio default intensity is a derived process in our model. It is shown to be a 
jump-diffusion that depends on the default level. This dependence is governed by a set of 
multiplicative factors - the so-called contagion factors. These factors enable a convenient 
semi-parametric specification which can be calibrated to any set of consistent portfolio 
tranche quotes. Furthermore, the fact that the model is two-dimensional and Markovian 
means that efficient lattice implementations are available. 

Our framework can be viewed as a low-dimensional short-rate version of the approach 
described by Schonbucher in [2]. In particular, Schonbucher describes a Heath- Jarrow- 
Morton (HJM)-like model of forward default transition rates which will in general be non- 
Markovian and require Monte-Carlo simulations. The BSLP model, on the other hand, 

1 Here the "top" refers to the portfolio-level dynamics, and the "down" refers to the single name 
dynamics which, if needed, might in principle be constructed consistently with the top-level portfolio 
dynamics. In this paper, we will only be concerned with the portfolio-level dynamics. 
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imposes local drift conditions that allow for fast lattice implementations, as indicated 
above. 

Both models coincide in the local intensity limit where the intensity process becomes 
a locally deterministic function of the default level. In this case, we obtain a continuous- 
time, one-dimensional Markov chain driven by a deterministic Markov generator. This 
can serve as a model in its own right, useful for constructing an arbitrage- free interpolation 
of liquid index (e.g. iTraxx or CDX) tranche prices across strikes and maturities. The 
interpolation can then be used to price non-standard index tranches consistently, which 
is often a problem in the standard base correlation framework. 

This is similar to the way local volatility models are used to price exotic equity or FX 
options off liquid vanillas. While such a framework can formally be viewed as a dynamic 
model, it is known that local volatility models misspecify the dynamics of the volatility 
process, and are therefore ill-suited for pricing instruments sensitive to the forward smile 
dynamics, such as forward starting options and cliquets. This points to the need for 
stochastic volatility models. 

Similarly, what is needed in the present context, for the pricing of multi-period credit 
derivatives is a stochastic evolution of the loss distribution. This is obtained by making 
the Markov chain generator stochastic. In this way, we arrive at the full-blown, stochastic 
intensity version of the BSLP model. 

Besides Schonbucher [2], our model is related to a few other models suggested pre- 
viously in the literature. In particular, it can be compared to the Markov "infectious 
default" model by Davis and Lo |12j^) In their approach the portfolio default intensity is 
piecewise deterministic, and follows a pure jump process that jumps upon defaults. In our 
case we have a stochastic, rather than piecewise deterministic, portfolio default intensity 
as a result of introducing a diffusion component in the dynamics of the driving factor. 
In addition, and perhaps more importantly, we set up a semi-parametric framework that 
enables accurate and fast calibration to market data. (See Appendix C for a further 
discussion of the relation between the two models.) 

Our model can also be compared to the time-changed birth model of portfolio loss 
dynamics of Ding et al. j9]. They consider a linear birth process that has a self-affecting 
property (controlled by a single parameter) and is therefore capable of modeling credit 
contagion of credit losses in the portfolio. To have more flexibility in the dynamics, an 
independent parametric time-change process (CIR or quadratic Gaussian) is introduced. 
This plays a very similar role to the driving intensity in BSLP. The main difference between 
the models is that BSLP is represented by a non-linear death process and has a semi- 
parametric specification with a level-dependent amount of contagion controlled by a set of 
contagion factors. Our semi-parametric approach enables accurate calibration to any set 
of consistent quotes. However, this comes at the price of losing the analytical tractability 
of Ding et al. and necessitating the use of numerical (or approximate analytical) methods 
instead. An additional difference between the two approaches is that, as indicated above, 
we admit a back-impact of defaults on the driving intensity process, which is absent in 
the model of Ding et al. 

2 which to our knowledge can be considered the first dynamic model of the top-down type, without 
the "down" part. See also Frey and Backhaus [1 3 a for a related approach to the portfolio dynamics with 
a mean- field interaction between individual defaults. 
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The paper is organized as follows. In Section [2] we present the ID local intensity version 
of the model. Section [3] describes the 2D stochastic intensity extension. We specify the 
default intensity process, forward equations for the full 2D model, and its calibration. 
In Section [4] we summarize lattice-based pricing algorithms for tranche, tranche options, 
forward-starting tranches etc. Section [5] contains detailed analysis of numerical results, 
including calibration of the model and pricing of index tranchlets and options on indices 
and tranches, as well as analysis of conditional forward spreads implied by the model. 
Furthermore, we explain how to calculate sensitivities in our framework, and compute 
index deltas for tranches and tranche options. Section [6] summarizes. In addition, four 
appendices discuss more technical model details. Appendix A derives the relation between 
the next-to-default intensity and the portfolio default intensity. Appendices B and C 
develop the continuous-time formulation of BSLP under different assumptions on the 
stochastic driver dynamics. Appendix D describes the adiabatic approximation to BSLP 
that can be used to construct a semi-analytical approach to pricing credit vanilla and 
exotic products. 



2 BSLP with Local Loss Intensity 

In this section we describe a local intensity version of the BSLP model, where default 
intensities are assumed to be locally deterministic (dependent on the loss level only). 
More specifically, we construct a Markov portfolio default process whose marginal default 
distributions will be consistent with any set of arbitrage-free quotes for tranches on the 
portfolio. Our construction is similar to the one-step Markov chain construction proposed 
by Schonbucher [2]. 

The model presented here will form the basis of the full-blown BSLP model, which 
will be introduced in Sect. [3] as a stochastic generalisation of the framework developed 
here. 

2.1 Markov chain for portfolio default process 

We construct a model for the default counting proces^] N t : 

N 

i=i 

where T{ is the default time of the i th name in the portfolio, and N is the initial number 
of portfolio names. The indicator function is denoted by I. 

We assume that the counting process Nt is Markovian, and model it as a continuous- 
time Markov chain with generator matrix A t . We recall that off-diagonal elements of 
At = { a ij(^)}f^=o provide transition probabilities for infinitesimal time intervals dt: 

P[N t+dt = j\N t = i] = aij (t)dt (2) 

for j 7^ i. 



3 Thc relation between Nt and the portfolio loss L t is described in section 2.2 
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Given the generators A t , the matrix of finite-time transition probabilities P(Nt, T\N t , t) 
with matrix elements 

Pij (t,T) = P[N T = j,T\N t = i,t] 
satisfies the forward equation 

PA (3) 



dP 
df 



In what follows we restrict our analysis to a piecewise time-homogeneous continuous- 
time Markov-chain, where different generator matrices are applied for time intervals 
[0, ti], [ti,t 2 ] etcQ Note that for any interval 1 2 ] o n which A t is constant, A t = A, 
the forward equation can be solved to give: 



Pij(ti,t 2 ) = (■ 



-h)A\ 

' ii 



(4) 



Without the piecewise-homogeneity of A t , we would need to use a time-ordered exponen- 
tial in the equation above, along with the substitution (t 2 — t\)A — > f t * A s ds. 

Clearly, default counting is a non-decreasing process, which means that its generator 
should have zeros below the diagonal Furthermore, assuming that at most one default 
can occur in the infinitesimal time dt, the generator can be taken to be bi-diagonal. Since 
we are dealing with infinitesimal time-steps this is a reasonable assumption, which makes 
a parsimonious description of the model possible^] 

The simplest Markov non-increasing process with a bi-diagonal generator is a well- 
known linear death proces^j The constant generator A takes the following form 



A = X 



( -N 




V 



N 
-(N-l) (JV-1) 













\ 



-1 1 

0/ 



(5) 



Here A is a parameter interpreted as the average single-name default intensity of all names. 

It is well know (see e.g. Feller [15]) that the death process gives rise to the binomial 
distribution for the number of events (in our case, defaults): 



7r(n, t) = P[N t = n}= (^) (1 - e~ xt ) n e 



-(N-n)Xt 



(6) 



The binomial distribution ^ corresponds to the case of non- interacting ( "uncorrelated" ) 
credit exposures in the portfolio, and is very far from distributions implied by market 
prices of liquid tranches, for any choice of A. 



4 Typically, h,t 2 , 
section 



2.21 



will be chosen to be maturities of the tranches that we are calibrating to (c.f. 



5 Equivalently, if instead of the number of defaults N t we model the number of surviving firms n t — 
N — N t , then the generator will have zeros above the diagonal. 
6 Multiple instantaneous defaults have been explored in [5]. 

7 Note that usually death processes are formulated in terms of the decreasing survival counter variable 
n t — N — Nt rather than the default counter variable Nt ■ 
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To produce a more interesting and non-trivial distribution of defaults, we consider a 
minimal modification of the homogeneous linear death process (|5]). In particular, we make 
it piecewise time-homogeneous to calibrate to tranche quotes at different maturities. We 
also generalize it to a non-linear death process, where transition intensities are a general 
function of the number of defaults n in the portfolio. We therefore consider the following 
generator: 

/ -fo(t)N f (t)N ••• \ 



A t = X 



\ 










-fi(t)(N 














-f N -l(t) 








fN-l(t) 





(7) 



/ 



Clearly, the linear death process ^ is recovered if all fi(t) = 1. The generator (m) 
corresponds to the case where the next-to-default intensity in the scenario where n names 
have defaulted is given by: 

\ N tD(n,t) = Xf n (t)(N-n) (8) 

If f n (t) 7^ 1, it means that the default intensity is a non-trivial function of the default 
level and the resulting implied default dynamics is contagious (self-affecting). Hence, in 
what follows, the parameters fi(t) introduced here will be referred to as contagion factors. 

As discussed in the next section, a suitable parametrization of the contagion factors 
will enable calibration to any set of consistent portfolio tranche prices. 

The bi-diagonal generator provides a convenient low-dimensional parametrization of 
the portfolio default process. Indeed, numerical calculations of break-even tranche spreads 
according to equations (11) and (14) (see below) require evaluation of loss distributions 
at quarterly time steps. Probabilities of multiple defaults within such periods cannot be 
neglected. Attempting to directly parameterize these probabilities in a discrete Markov 
chain framework would be a difficult task. 

In contrast, the continuous-time framework allows the calculation of all finite time, 
multiple default, probabilities by taking matrix exponentials of the generator ([7]). This 
can be done efficiently if we neglect small differences in accrual periods, as described in 
the next section. 



2.2 Numerical implementation and calibration to tranche quotes 

We assume that we are given quotes for a set of tranches at different maturities. Let 
k denote the ordered set of all low strikes (including 100% if we are calibrating to the 
index level). For example, for iTraxx we have k = {0%, 3%, 6%, 9%, 12%, 22%, 100%}. 
Similarly, the set of quoted maturities is denoted by T, so that for example for iTraxx 
series 6 we would have: f = {20-Dec-09, 20-Dec-ll, 20-Dec-13, 20-Dec-16}. 

Since tranche payoffs are determined by the cumulative portfolio loss, L t , and not the 
number of defaults, N t , we need to specify the relation between the two variables. The 
most general assumption we can make within our model is that L t is a function of N t 
onl)|^J i.e. L t = L(N t ). Typically, however, we can set the recovery rate, R, to a constant 

8 In particular, this excludes time-dependent recovery rates. 
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value, in which case we have: L(N t ) = (1 — R)N t /N, where N is the initial number of 
names in the portfolio. 

In order to calibrate to the given set of quotes we need to choose a parametrization of 
the contagion factors fi(t). To do this we introduce a contagion function, g(L t ,t), which 
is a function of the portfolio loss level and time. This function is defined on the interval 
[0, 1] x [0, oo) and takes values on the positive real line: [0, oo). We define the contagion 
factors as: 

fi(t) = g(L(i),t) (9) 

where as above L t = L(N t ). Lastly, we restrict the functional form of g so that the 
number of free parameters exactly matches the number of tranche quotes. To do this we 
assume that g(l, t) is piecewise constant in t (for a given /) with node points given by the 
quoted maturities, T. In addition we take g(l, t) to be piecewise linear in I (for a given t) 
between node points k (and piecewise flat outside). 

Using a multi-dimensional solver we can now calibrate the contagion factors /&(£) to 
the observed tranche prices. This is done by calculating all finite time loss distributions 
and using them along with the standard pricing formulae (reviewed in the next section). 
This results in a calibrated generator of the default process: 

a l3 (t) = X(N - i)f i (t)(5 i+hj - S itj ) i, j = 0, . . . , N (10) 

where 8^ denotes the Kronecker delta. As described previously, finite time transition 
probabilities can be obtained from the generator ^ by matrix exponentiation. This can 
be done using efficient numerical methods such as the Pade approximationR 



2.3 Pricing 

Once all finite-time loss transition probabilities are calculated from the generator, the 
pricing of a tranche with lower strike Kd and upper strike K u (both defined as percentages 
of the total portfolio notional Nq) with maturity T is done using standard formulae which 
we present here for the sake of completeness. 

In the following, we assume that the recovery rate R, is fixed for all names in the 
portfolio. Then the default (or contingent or protection) leg, £>, of a tranche paid by the 
protection seller to the protection buyer is given by 




where the sum runs over all coupon dates Tj, i — 1, . . . , M, B{ = 5(0, T^) is a risk-free 
discount factor, and ELi is the tranche expected loss, given by: 

5L^EL T! =E[L [M] (T,)] (12) 

9 In practice, we can neglect small differences in accrual periods, and use the Pade approximation to 
calculate the transition matrix Pi for the first period only. All other transition probability (two-period 
etc.) matrices are obtained by repeated multiplication of P\ by itself. 
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where L[K d ,K u ](t) is the tranche loss at time t given by: 



L [KM (t) = {L t - K d )+ - (L t - K u ) + (13) 
The premium leg (paid by the protection buyer to the protection seller) is given by 

M / p t- 



V(S) = S-N Y,^i[B i -EN Ti - J t T —^- B(0,t)dEN t ) 

M 

~ S-N J2^B t -(EN t ^ + EN t ) 



(14) 



8=1 



where S is the tranche spread, Aj is the day count fraction and ENi is the expected 
tranche outstanding notional at time Tj. It is defined by: 

ENi = EN Ti = E[N [Ku , Kd] (Tj] = E[K U -K d - L [KM ] (15) 



The integral term in (14) represents the accrued coupon due to defaults happening 
between the coupon payments dates. The integral is calculated using the standard ap- 
proximation: (t — Tj_i)/(Tj — Tj_i) — > 1/2 and B(0, t) — > 

The tranche par (or fair or break-even) spread, S, is determined from the par equation 
T> = V(S). Most index quotes are given in terms of the par spread. For the equity tranche, 
however, the market convention is to charge an upfront payment from the protection buyer 
while fixing the running spread S at 500 bp. 



2.4 Bootstrapping 

It can be useful in practice to convert the calibration into a bootstrapping problem. This 
is possible if we use piecewise-flat interpolation of the contagion function g(l,t) in both 
time and loss dimensions. 

This can be easily seen by the fact that the calculation of the par-spread, S[K d ,K u ] {T) 
(of a tranche with low strike K d , high strike K u , and maturity T) will depend only on the 
function g(l, t) in the region < t < T and < / < K u . 

Because of the piecewise-flat assumption we can calibrate the contagion function values 
at the points k and T iteratively in the time and loss dimensions. 

At this point we should note why we have not made the piecewise-flat loss interpolation 
choice in the previous section. This is because for typical tranche quotes such as iTraxx, 
the piecewise-flat assumption is not very good for the equity region and results in poor 
calibration. This is essentially due to the thickness of the standard tranches, e.g. 3% for 
junior iTraxx and CDX tranches. Indeed, if we calibrate to thinner tranches then the 
piecewise-flat assumption will become more viable. 

In particular, in the limit where we consider tranches that are exactly one default wide, 
piecewise-linear and piecewise- const ant interpolations become equivalent. This means 
that bootstrapping will produce very accurate calibrations provided the input quotes are 
arbitrage free. 
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2.5 Discussion 



In this section we discuss the applicability of the local intensity model presented above. 

First we note that the price of a tranche depends only on the marginal default distri- 
butions at a set of time points but not on the inter-temporal correlations between defaults 
at different times. In other words, tranches are single-period instruments. 

Let us hence assume that we have found a set of marginal default distributions such 
that prices of all liquid tranches are matched. It is well known that this is not equivalent 
to, and not by itself sufficient for, defining a dynamic model. 

This is best seen by considering the path probability for the default counting process 
N t to pass through levels [N h , N 2h , • • • , N nh ] at times t = [h, 2h, . . . , nh]. Using the chain 
rule, we write 

P [N h , N 2h , ...,N nh ]=P [N h ] P [N 2h \N h ] ...P [N nh \N h N 2h . . . N {n _ 1)h ] (16) 
The important point here is that tranche market prices do not constrain the conditional 



transition probabilities in (16). The only thing that can be uniquely determined from 



the marketpj are marginal loss distributions at various times (implied distributions) and 
one-step transition probabilities 11 P[Nt+h\Nt\. 



The calculation of the path probability (16) reduces to a product of one-step transition 



probabilities only in the special case where Nt is a Markov process. However, whenever 
additional stochastic risk factors are present in the model this would generally not be the 
Hence it is only by making further assumptions that a set of marginal distributions 



case 



12 



for different time horizons, originally viewed as a disjoint set, can be "glued" together 
within a unique dynamic process. 

In the local intensity setting, we assumed that the portfolio default counting variable, 
Nt, is the only relevant state variable, and its dynamics was further supposed to be 
Markov. This assumption allowed us to construct a model which can calibrate perfectly 
to any set of tranche prices. 

Formally, such a ID Markov chain model may be viewed as a dynamic, multi-period, 
model since the Markov property allows one to calculate all default transition probabil- 
ities. However, the assumptions underlying this model are likely to be too extreme. In 
particular, defaults are rare events and we expect portfolio spreads to move stochastically 
even in the absence of defaults, which cannot be achieved in the present ID setting. 

Even if we are willing to ignore the requirement of any realism in modeling portfolio 
spreads, we face the problem that the default-only model is completely specified once 
calibrated. It is likely that we will need extra flexibility to adjust the dynamic properties 
of the model when pricing exotic derivatives. 

For example, stochastic loss intensities are needed if we want to have control over 
the tranche option implied volatilities. Another issue is that the ID model implies high 



10 Assuming its completeness, i.e. the availability of tranche price quotes for all strikes and maturities. 

11 This was demonstrated by Britten- Jones and Neuberger (BJN) [T3] in the context of stochastic 
volatility models for equity options. Modifications of their derivation for the present context of tranche 
pricing are straightforward. 

12 Let M t = (Nt,...) be the vector combining all relevant state variables. If the multi-dimensional 
process M t is jointly Markov, then its individual components N t . . . are generally non-Markov when 
viewed separately. 
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inter-temporal correlations between losses at different time horizons. Again we need to 
introduce stochastic intensities to gain more control. 

Furthermore, experience with equity and FX modeling suggests that local volatility- 
type models tend to give rise to strong time-inhomogeneity of dynamics which, while 
allowing for a good fit to today's option prices, makes the evolution of option prices (or 
equivalently, implied volatility) unrealistic. 

For these reasons the ID model is unlikely to be adequate for pricing path-dependent 
instruments or/and derivatives that are sensitive to the dynamics of the loss process. Note 
that most exotic credit derivatives (tranche options, forward- starting tranches, leveraged 
super-seniors etc.) are of this type. To model these instruments, we need a dynamic 
extension of the framework developed so far. This will be presented in the next section. 

To summarize, the local intensity version of our model appears inadequate for pricing 
credit exotics, however it is useful on two other counts: 

• It provides a very convenient way to parameterize a set of observable liquid tranche 
quotes. This can be used for an intuitive interpolation-based pricing of non-standard 
tranches off the liquid ones. 

• It can serve as a first step towards the second, dynamic stage of the model. 



3 Dynamic Intensity BSLP model 

3.1 Stochastic loss intensity 

We construct a dynamic version of the model presented above by promoting the constant 
A in the generator dTl) to a stochastic driving process Y t . This means that the next-to- 
default intensity \mD(Nt,Y t ,t) is now stochastic and is given by: 

Xm D (N t = k,Y t ,t) = (N- k)F(k,t)Y t (17) 

Here we have introduced the contagion factors F(k, t), which, as before, are a deterministic 
function of default level k and time t. Note that these factors will, in general, differ 
numerically from the factors / defined in equation (10). We will occasionally refer to the 
former and the latter as the 2D and ID contagion factors, respectively. 
We consider the general form of the SDE of the driving process Y t : 

dY t = fjt(y t , t)dt + a(Y t , t)dW t + 7 (y t , t)dN t (18) 

Several specifications of the dynamics of Y t are possible. We have considered a mean- 
reverting, log-normal process with 

tii n (Yt-,t) = a y t _(ln0 t -]ny t _) (19) 
ai n {Y t _,t) = oY t _ 

lln{Y t _,t) = 7 Ft_ 
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and a mapped OU form Y t = q(yt) where q(x) is a function ensuring positivity of the 
mapping and yt is an Ornstein-Uhlenbeck (OU) process with 



fJ>ou{yt-,t) = a(9 t -y t -) 

a ou (yt-,t) = a (20) 
7ou(yt-,t) = 7 

It will be important later that the Y t process does not depend on the default level 
N t . It can, however, include a jump component, dN t , which depends on changes in the 
default levefB 



Note that for both (19) and (20) the model obtained in the limit a — > 0, a — > oo and 
7^0 becomes equivalent to the local intensity model described in the previous section. 
Alternatively, in the zero volatility limit o — > with 7 ^ 0, we obtain a pure jump 
process for Y t driven by jumps in N t , which fits into the piecewise-deterministic Markov 
framework of Davis and Lo (DL) 11 1 



3.2 Index swap intensity process in BSLP 

We now consider in more detail the BSLP dynamics implied by our definition of the next- 



to-default intensity (17). This is done most conveniently by considering the SDE for the 
index swap intensity ^ X s (t). The reason for focussing on the swap intensity is that it 
represents the average portfolio single name intensity. In the absence of contagion this 
can be expected to be largely independent of the portfolio size. This is not the case for 
the next-to-default intensity, which even in the independence case is proportional to the 
remaining number of names in the portfolio. 

In what follows we assume that recovery rates are constant and homogeneous across 
the portfolio. If in addition we assume that defaults cannot happen simultaneously, then 
X s is obtained from the first-to-default intensity by division by the number of surviving 
name£3 

X s (t) = XmD N { ^ t) =Y t F(N t ,t) (21) 

While the relation between X s and XmD m ay look intuitively obvious, we present a formal 
derivation in Appendix A. 



13 Note that a slightly more sophisticated form of the jump term in ( 18 1 7/„(Yt, t) — jptY t or jou(Yt,t) = 
"fPt (where pt is a Bernoulli random variable) could be considered as well. We do not pursue this 
modification in what follows. 

14 The difference being the different number of possible states for Y t - in the original DL model, the 
loss intensity can only assume two values, while here it would be continuous until of course we discretize 
it for practical implementation. 

15 This is the stochastic intensity that prices the index correctly, assuming it is a single name CDS. For 
a homogeneous portfolio this is equal to the individual single-name portfolio intensities. 



16 Note that the relation (21) expressing the portfolio loss intensity X s (t) as a product of the driving 
factor Yt and the contagion factors F(Nk,t) resembles the relation between the stochastic time change 
rate v t and the generator of the initial Markov chain for the default process in the model by Ding et al 
[9], and can be interpreted correspondingly. 
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To derive the SDE we first note that using Ito's lemma for a pure jump process N t , 
we obtain 



dF(N t ,t) 



dF t (N t _,t) 

at 



dt + (F(N t ,t)-F(N t „,t))dN t 



(22) 



= F'(N t -,t)dt + AF(N t , t)dN t 



where N t _ = N t — 1 if there is a default at time t. Now the dynamics of X s (t) are 
determined as follows: 

d\ s (t) = d(Y t F(N t ,t)) = F(N t _,t)dY t + Y t .dF(N t ,t)+ 1 (Y t _,t)AF(Nt,t)dN t (23) 



Then (23) together with (22) and (18) yield the following SDE for the swap intensity 
A,(/): 

d\ a (t) = fi s (X s , N t , t)dt + <t.(\„ N t , t)dW t + 7 S (A S , N t , t)dN t 



(24) 



where 



fi s (X s ,N t ,t) 
(J s (X s ,N t ,t) 
j s (X s ,N t ,t) 



F{N t _,t)ii 
F(N t -,t)a 



A, 



F(N t -,ty 

A, 



t )+x. 



F'(N t -,t) 
: F(N t -,t) 



F(N t -,ty 

(F(N t -,t) + AF(N t ,t)) 1 



(25) 



A., 



F(N t .,ty 



t)+x, 



AF(N t ,t) 
F(N t -,t) 



Equations (25) determine the N t - dependent (and generally non-linear) transformation 
needed to obtain the coefficient functions of the SDE (24) from the coefficient functions 
of the SDE (18). It is interesting to consider these relations for the example of the log- 
normal mean reverting (Black-Derman-Toy or BDT) specification (19). For this case, we 
obtain 



fx s (X s ,N t ,t)/X s (t-) 
a s (X s ,N t ,t)/X s (t-) 
ls (X s ,N t ,t)/X s (t-) 



a (\n(F(N t -,t)9t)) - In X s (t-)) + 



F'(N t -,t) 
F(N t .,t) 



a 



(26) 



AF(N t ,t) 
7+ F(N t .,t) { +1) 



This is an interesting result. It shows that if we start with a jump-diffusion BDT speci- 
fication for the driving 1^-process, then the SDE for the quasi-observable swap intensity 
process X s {t) is obtained from that of Y t by adjusting its drift and jump intensity, while 
keeping the volatility constant. 

We can also clearly see that there are two qualitatively different types of contagion 
implied by the model. The coefficient 7 S (A S , N t ,t) determines the jump size of the swap 
intensity contingent on portfolio defaults. The jumps in intensity will mean-revert over 
time, hence we say that 7 is responsible for transient contagion. Note that even if we start 
with the pure diffusion case 7 = for the process Y t , the intensity X s (t) acquires a jump 
component with the jump size determined by the relative change of the 2D contagion 
factor F(N t ,t). 
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The second important effect of portfolio defaults is that they change the mean reversion 
level of the swap process. This is referred to as permanent contagion and is driven by the 
contagion factors F(N t ,t). 

This can be compared with the formalism suggested by Lopatin and Misirpashaev 
(LM) [TT]. In their formulation, they start with a stochastic next-to-default intensity 
\ LM (N t ,t) where the A^-dependence arises entirely through the drift term. Furthermore, 
LM do not admit jumps in the SDE for X LM (N t ,t), in contrast to the BSLP formula- 
tion. Hence, in the LM model next-to-default intensities do not jump upon defaults but 
gradually diffuse to a new mean-reversion level. 

However, jumps are present in the SDE for the swap intensity, which is related to the 



next-to-default intensity by the general formula (21): 



A^(0 = (27) 

Hence, we have: 

d\ L s M (t) = w L w to LU (N t ,t) + X^Nt^dj^j^ (28) 

To evaluate the second term, we use Ito's lemma for jump processes, assuming in addition 
that dN t <C iV — N t so that we can use an approximation: 

1 ( 1 1 \ .„ / 1 x " 



d = )dN t ~[ dN t (29) 

N — N t \N-N t N-N-J \N-N t J v ' 

This means that in the LM model jumps in the portfolio swap intensity \ s (t) are solely 
due to the reduction of the portfolio as defaults occur. In addition the jump size scales 
down in a fixed manner as (N — N t )~ 2 . This is in contrast to the BSLP model where the 
jump size of the swap intensity depends on the amount of contagion in the data, expressed 
via the 2D contagion factor F(N t ,t). 



3.3 Forward equation and transition probabilities in 2D 

We postulate that the 2D dynamics of the pair (Nt, Y t ) in BSLP model is Markovian. We 
find it convenient to use a simple discrete-time formulation of the model in what follows, 
with a time step dt considered as small but finite. A continuous-time formulation of BSLP 
may still be of interest from both theoretical and practical viewpoints. This is discussed 
in some detail in appendices B to D. 

Discretizing the range of Y t to a finite set {Y n }, the system is described in terms of 
the joint marginal probabilities 

n(j,n,t) = P[N t =j,Y t = Y n ] (30) 

and conditional transition probabilities 

p jm \kn(t, t + dt)=P [N t+d t = k, Y t+dt = Y n \N t = j, Y t = Y m ] (31) 
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The forward equation takes the form 



(k,n,t + dt) = ^Pjmiknjt, t + dt)ir(j, m, t) (32) 



7T ^ 



We note that due to the composition law of probabilities we have the following relation 

Pjm\kn(t,t + dt) (33) 

= P [N t+dt = k\N t = j, Y t = Y m ] P [Y t+dt = Y n \Y t = Y m , N t = j, N t+dt = k] 
The form of the second factor 

P mn (j, k,t)=P [Y t+dt = Y n \Y t = Y m , N t = j, N t+dt = k] (34) 



is fixed by the SDE ([18]). Note that here we used the fact that our SDE for Y t does not 
depend on the default level N t but only on the change in default levels dN t . 

The formulae presented so far are completely general. BSLP is defined by the func- 



tional form of the first factor in (33). 



g jk (Y m , t)=P [N t+dt = k\N t = j, Y t = Y m ] (35) 



This is determined by our definition of the next-to-default intensity in equation (17) as 
follows: 

( dt(N - N t )Y m F 3 {t) + 0{dt 2 ) ilk = j + l 
g jk (Y m , t) = l 1 - dt(N - N t )Y m F 3 (t) + 0(dt 2 ) Hk = j (36) 
y 0(dt 2 ) otherwise 

where we have introduced Fj(t) = F(N t = j,t). 
3.4 Tree discretisation 

To turn the above into a practical scheme, we discretise BSLP on a tree. To do this we 
introduce a discrete timeline to, t\, . . . , t n with (finite) time step At = £ i+1 — U. 
We assume that Y t and F^(t) are piecewise constant between timeline points. 



The transition probabilities defined in equation (36) are given to first order in dt. If 
we wanted to use these equations directly this would require very small (daily) time steps 
in the tree. 

In practice our discretisation step At is likely to be larger. In this case we can no longer 
assume that there is only one default over the time period. Our conditional transition 
probabilities are now given by: 

P[N t+ At = k\N t = j, Y t = Y m ] = ( e AtY ™ A A (37) 

\ / jk 

where A t is defined as in (M) with the substitution A — > 1 and fi —> Fi(t). This expression 



generalizes (36) to the case of a finite time step At. Note that for a small but finite At, 



(37) coincides, to the first order in At, with (36) integrated over the interval [t, t + At] 
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assuming that Y t and Fk(t) do not change in this interval. In other words, (37) can be 
interpreted as the transition probability matrix in a conditional discrete-time Markov 
chain obtained by "freezing" the random variables Y s and Fk(s) for all t < s < t + At at 
their values at time t. 



Given equation (37) we can construct a discretisation of BSLP on a 2D tree. In partic- 



ular, the factorisation in equation (33) means that to calculate 2D transition probabilities 
we can proceed in two steps: 

1. Discretise the Y t process on a lattice or tree and calculate all ID y-transition prob- 
abilities. 



2. Calculate the conditional transition probabilities in (37). As mentioned previously, 
efficient numerical methods, such as the Pade expansion, for calculating matrix 
exponentials are available. 

The probabilities Pj m \kn(t,t + dt) are then given by the product of the probabilities cal- 
culated in the two steps, as in (33). 



3.5 Two dimensional calibration 

BSLP is specified by the choice of contagion factors, F(N t , t) and the parameters governing 



the Yj-process (18). As such there are two possible routes for calibrating the model to 
external data. 

We could fix the form of the 2D contagion factors F(N t ,t), and try to fit the Y t to 
liquid tranche prices. Alternatively, we could postulate some fixed law for Y t and calibrate 



the model in terms of the F(N t ,t). In this paper, we explore the second route 17 

In more detail, once we have fixed the parameters that govern the y r process, we can 
calibrate the contagion factors F(N t ,t) exactly as we do in Sect. [2j In particular, we 
can still use the bootstrapping algorithm as described. The only difference is that now 
tranche prices are calculated on a 2D tree. This makes calibration slower, although the 
speed is still acceptable (about 1-2 min for a standard index portfolio with calibration to 
four standard maturities). 

However, a much faster calibration algorithm is also available. This is described in 
the next sections. 



3.6 Fast calibration by forward induction 

In this section we present a fast calibration algorithm that enables a re-calibration of the 
full 2D model starting from a calibrated ID default tree (lattice) constructed in section [2j 
It uses a recursive procedure of "integrating in" the stochastic loss intensity. Our method 
is similar to the algorithm developed by Britten- Jones and Neuberger (B JN) in the context 

17 Note that Lopatin and Misirpashaev [TT] essentially choose the first route by making their analog 
of the Yf-process explicitly dependent on the A^-variable, while effectively taking the factors F(N t ,t) 
to unity. Computationally, the two routes turn out to be nearly identical, however as we have seen the 
difference in parametrization does matter as it leads to observable effects in the dynamics of the swap 
intensity. 
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of stochastic volatility modeling. In turn, the B JN method is closely related to the Markov 



projection method used by Lopatin and Misirpashaev pj. 18 

We start by summing over n in the forward equation (32), to get the probability 
unconditional on Y t+dt : 



ir(k, t + dt) = yP [N t+dt = k\N t = j, Y t = Y m \ ir(j, m, t) 



(38) 



This is the marginal iV-probability at time t + dt. Assuming that the model is constructed 
as described above (i.e. in two stages) we have the ID iV r tree (lattice) calibrated to the 
observed set of tranche quotes. Marginal probabilities in the full 2D model should match 
ones calculated in the ID model. Moreover, we can relate the marginal probability at 
t + dt with the one at t using the forward equation in the ID model: 



%(k,t + dt) = ak-i,k^(k — l,t)dt + (1 — ak,kdt)n(k,t) + 0(dt 2 ) 



(39) 



We now re- write equations (36) in the following form: 



P[N t 



t+dt 



k\N t =j,Y t = Y rr , 



dt(N - 3)Y m q ] {t)f, J {t) + 0(dt 2 ) 
1 - dt(N - ] )Y m q j {t)f, J {t) + 0(dt 2 ) 
0{dt 2 ) 



if k = j + 1 

if k = j (40) 

otherwise 



where fj(t) = 
Qj(t) = q(N t 



f(N t = j, i) are the contagion factors introduced in equation (10), and 
= j,t) are the drift adjustment factors. 



We assume that the contagion 
factors fj(t) are chosen such that the ID iV r tree is calibrated to the market data. Note 
that as long as the factors qj{t) are not yet defined, the parametrization in (40) is only a 
matter of convenience, and is completely equivalent to (36). 



By substituting (40) into (38) and (10) into (39) and by comparing the two resulting 



forward equations, we obtain the following constraint on the drift adjustment factors qj(t): 



A 



E[y*|iv t ; 



(41) 



The no-arbitrage drift constraint just derived is our short-rate analog of the no- arbitrage 
drift constraints in the H JM-like construction of Schonbucher [2] , where the drift correc- 



tions are typically non-local in time. In contrast, our drift constraints (41) are local in 



both N t and t, and are thus amenable to a fast lattice implementation. Note that our 
drift constraint coincides with the dt — > limit of the drift constraint obtained by BJN 
[H] using a more involved argument (and in a different context). Also note that the 
same drift condition (albeit with a different notation) is obtained by LM [11] who further 
identify it as a special case of a general Markovian projection approach. 



18 When the BSLP model was developed, we were not aware of V. Piterbarg's paper on the Markov 
projection method (V. Piterbarg, "Markovian projection method for volatility calibration" , available at 



http://ssrn.com/abstract| = 906473), and instead referred to this approach as the BJN method. It appears 
that, at least for the setting used in this paper, the two methods are very similar, if not identical. Note 
that while Piterbarg does not cite BJN, both BJN and Piterbarg cite the work by Dupire on the link 
between stochastic and local volatility models. Dupire 's approach seems to provide a common basis for 
both the BJN and Markov projection methods. 
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We now use the BSLP drift constraint (41) in order to set up a convenient and fast 



calibration method for the 2D BSLP tree based on a combination of the ID calibration 
of the iVt-tree and a forward induction method. 

We assume that the ID calibration of the iVVtree is performed as discussed above. We 
start with the initial conditions for the 2D and ID probability distributions, correspond- 
ingly, 

7l(N , Y m , 0) = 5v ,o5m,m , ^{N , 0) = 6 Nofi (42) 
where rh is the index corresponding to the initial value Yq (which we assume to be known). 



Using Eq.(41), we solve for q No (0): 



A7r(iVo,0) A , , 

^o(O) = v Y y N } y m = tt , No = (43) 

Note that for No ^ the correction factors gv o (0) are undefined. However, this does not 
pose any problem as such states are unachievable at time t = 0, and therefore play no role 
in the dynamics whatsoever. If desired, these parameters can be assigned some definite 
dummy values that would not have any impact on the numerical results reported below. 
We now use the forward equation 

n(k,n,t + dt) = ^ 9k-ijs(Y m , t)P mn (k - 1, k, t)ir(j, m, t) (44) 

m 

+ ^2gk,k( Y rrn t ) P mn(k,k,t)7l(k,m,t) 



(where the conditional Y t - transition matrix P mn (i,j, t) is defined in (34)) at t = to calcu- 



late the joint probability distribution ir(k, n, dt). Then we calculate the drift adjustments 



qj(dt) for the second period from the drift constraint (41). Next we use it to calculate 



Tv(k, n, 2dt) from the forward equation (44) evaluated at t = dt, and so on. As a result, we 
have a full 2D BSLP tree calibrated to the set of tranche quotes using a fast and effective 
algorithm. 

To avoid possible misunderstanding, a comment on the meaning of the above procedure 



is in order here. The drift condition (41) provides a way to specify the 2D dynamics 
once the ID dynamics of N t is fixed. Is this equivalent to saying that the 2D dynamics is 
uniquely (once the SDE for Y t is specified) fixed by observed tranche prices? The answer is 
no, of course, as market incompleteness means there is no unique correspondence between 
tranche prices and marginal default distributions. The meaning of the above procedure 
is rather to pick a ID model that matches all tranche quotes, and then calibrate the full 
2D dynamics to this model, not to the data directly. 



3.6.1 Fast calibration in practice 

The algorithm of forward induction-based calibration of the BSLP model just presented 
is simple and intuitive, however it is not ideal from a practical viewpoint, as it assumes 
that the time steps are small enough to justify the use of a binomial approximation for 
the one-step default process. In practice, this means we should take daily steps, which 
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may slow down calibration and pricing. If we want to be able to have a tree with larger 
steps (e.g. monthly), we need to use a different method. 

Several alternatives of different complexity can be considered at this point. The first 



one is to generalize the relation (40) to the case of non-infinitesimal time steps At. We 



have considered the following specification 

Q jk (t) {e MY - M ) ]k Hk>j 

P [N t+At = k\N t = j, Y t = Y m ] = { 1 - P [N t+At = k\N t = j, Y t = Y m \ if k = j (45) 

otherwise 



which is similar to parametrizations used for stochastic volatility models by BJN 
Here A t stands for the calibrated generator of the ID default-only model, and Qjk{t) 



are drift adjustment factors similar to the factors qj(t) in (36). We may interpret the 
term e AtYmAt in (45) as an "initial guess" or a "prior" for the finite-time conditional 



transition probability, which is then corrected by a set of multiplicative drift adjustment 



factors Qjk using a finite-time version of (41). Note that the ansatz (45) is somewhat 
non-symmetric with respect to its dependence on Y m , i.e. the drift adjustments Qjk 
with j ^ k are assumed to be independent of Y m , while the diagonal adjustment Qjj is 
implicitly dependent on Y m , as required by conservation of probability. 

We have implemented this scheme and tested it on a number of portfolios. Unfortu- 
nately, we have found that while this method works well for some portfolios, it develops 
numerical instabilities for others. However, we were able to find a simple practical solu- 



tion to this problem, which involves using the adjustment given in the first line of (45) 
for both off-diagonal and diagonal transitions, and then rescaling all probabilities by a 
common factor chosen to ensure the correct normalization. Note that this produces a 



more uniform dependence on Y m than in the ansatz (45). This algorithm was found to be 



stable and accurate in all cases we testec 19 , with model outputs similar to those obtained 
with a direct 2D calibration described above. 

Lastly, we note that in addition to the numerical methods discussed so far, analyti- 
cal approximations to the model are possible as well. In particular, we can consider the 
adiabatic approximation, which is expected to produce accurate results as long as the 
characteristic time of changes in spreads is much smaller than those in the loss counting 
variable. This assumption is expected to hold in reality (spreads change daily, while de- 
faults are rare events). Derivation of the adiabatic approximation is presented in Appendix 
C. Interestingly, this approach produces an analytical formula for the drift adjustments 



Qjif) similar to the one defined by (41) 



4 Pricing algorithms 

Models discretised on a tree or lattice are particularly suitable for pricing products which 
have a payoff that is amenable to a backward recursion algorithm. Typically, these are 
products with embedded optionality, such as tranche options or leveraged super-seniors. 

19 The price we have to pay with this method is that it introduces a small mismatch between tranche 
prices evaluated in the 2D and ID versions of the model, but mismatches were found to be negligibly 
small for all test portfolios we considered. 
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In addition, tree pricing^ of products that are weakly path-dependent is also tractable. 
By weakly path-dependent we mean that the payoff depends on the loss path at only a 
few points in the past. A forward tranche is an example of such a product. 

In this section we take a closer look at how tree pricing algorithms can be applied to 
the products mentioned above. Note that the algorithms presented are not specific to the 
BSLP model. In the following we also assume that defaults and losses are related simply 
by Lt — (1 — R)Nt and hence we will only talk about losses below. 



4.1 Tranche pricing by backward induction 

Let T>(t) and V(t) be the default and the premium legs at time t of a tranche with 
strikes Kd and K u with maturity T. Let = t%, . . . ,tM = T be the time grid and T n 
(n = 1,2, ... , M c ) be the coupon payment dates on the grid. Then, neglecting the accrued 
coupon, we obtain the following expressions: 



E; 



E; 



M-l 



B{ti,t j+1 ) (L[ Kd)Ku ](t j+1 ) - L[ KdtKu ](tj)) 

j=i 

K d ,K u ] 



(46) 



where the tranche loss, L[ Kdi Kj, and tranche notional, N[ Kd)K i, have been defined in 
section 2.3 Here A n is the day count fraction for the period [T n _i,T n ] and Ej stands 



for the expectation value conditional on the history up to time ij. Assuming that the 
last time grid point falls on a coupon payment date, we impose the following boundary 
conditions for T> and V (here % is the A-index and n in a Y- index): 





^■MN [KdtKu] (t M ) 



(47) 



Using the tower law for conditional expectations 

E t [X T ] =E[X t \ H = E [E [X T \ H I T t \ 



t< s <T 



(48 



and splitting off the first terms of the sums in (46), we obtain the recursive formulas 

Lt 



V(U) = B(U, [V(t i+1 )] + B(ti, U+i) (Ei [L [KM (ti+i 

V(U) = S(* i ,t i+1 )E i [7'(* i+ i)]+A n JV [jr(( ,j ftl] (T ft ) < y tjl r B 



J [K d ,K % 



(49) 



Here the second term in the second equation is non-zero only on coupon payment dates, 
and corresponds to a coupon added on that date. Using the one-step backward equations 
to evaluate the expectations entering (49), we obtain the backward induction pricing 
method for a tranche. 



20 It is also possible to construct a Monte Carlo implementation of BSLP to price path dependent 
products. This is not explored in this paper. 
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4.2 Tranche option 



Given the algorithm for calculating the tranche mark-to-market by backward recursion 
on a tree, we can price a tranche option using the standard tree option pricing procedure. 
Note that today's value, Vq of the tranche option with strike k and exercise date T\ and 
maturity T 2 is given by: 

V = E [M Tl (T2,k) + } (50) 

where A4t(T, c) is the mark-to-market at time t of the underlying tranche paying coupon 
c and with maturity T. Therefore, to price the tranche option, we build a tree or lattice 
up to maturity T 2 , and then roll the tranche price backward on this tree from T 2 to T\. 
This provides the boundary condition at T\ for the tranche option which is then rolled 
backward in time from t = T\ to t — 0. 



4.3 Forward starting tranche 

Forward starting tranches provide protection against tranche losses in a pre-specified 
future period [t,T]. The distinguishing feature is that defaults occurring prior to t do 
not affect the subordination of the tranche ^_ In particular, a forward tranche with low 
strike Kd and high strike K u can be valued as the forward value of a tranche with adjusted 
strikes, K' d and K' u : 

V = E [M t (T;K' d7 K' u )} (51) 

where Ait(T; K u , Kd) is the mark-to-market at time t of a tranche with maturity T and 
low and high strikes Kd and K u . The strikes are adjusted by the loss, L t , at time t and 
given by: K' u = min(l, K u + L t ) and K' d = min(l, Kd + L t ). This dependence of the payoff 
on the loss makes the forward tranche path dependent. 

To price this on a tree we need to know the mark-to-market, A4i(t), of the tranche at 
all loss nodes, i, of the tree at time t. The payoff of the forward tranche is dependent on 
the loss at t and hence each A4%(t) needs to be calculated on a separate sub-tree emanating 
from node i. The values A4i(t) provide a boundary condition for the tree at time t which 
can be rolled back to t = as usua0 



4.4 Leveraged super senior 

Leveraged Super Senior (LSS) is a tranche with the added feature that the trade knocks 
out once a certain trigger is breached. Different versions define the trigger to be either 
the portfolio loss or spread, or MTM of the tranche. In particular, for the loss trigger, the 
trade knocks out when the portfolio loss L t hits a (deterministic) barrier L(t) for the first 
time. Assuming a deterministic and fixed recovery R as before, we can translate this into 

21 Hence the forward tranche is not simply the difference of two standard tranches with maturities given 
by t and T (the latter is sometimes referred to as the plain vanilla forward tranche). 

22 This pricing algorithm can be somewhat simplified if the forward induction-based calibration of 
section [3TB] is used. In this case, marginal default probabilities are calculated at the calibration stage. 
Hence, we only need to roll the tranche MTM backward in time from T to t. The price at time t = is 
then given by a weighted sum of A4i(t) at all nodes i, with the weights being the state probabilities at 
time t. Note that the same comment applies to tranche options as well. 
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an equivalent default count boundary N(t) = L(t)/(1 — R). The random hitting time r 
is therefore the following: 

t = inf{t : N t > N(t)} (52) 
The payoff for the protection buyer is given by: 

P(t) = mm (K, M(t + At)) = M(t + At) - (M(t + At) - K) + (53) 

where K is the collateral posted and At is the unwinding period (typically two weeks) 
during which the trade is terminated and unwound after first breaching the barrier. Equa- 



tion (53) implies that the LSS can be viewed as a portfolio of a long position in the 
super-senior tranche and a short position in an American barrier option (the "gap risk 
option"). Let C(t) be the price of this option at time t. Assuming for simplicity that the 
length of the unwinding period At is equal to the time step on the tree, the option can 
be priced using the standard backward recursion: 

C{ti) = B(ti, ti+iJIf^jv-^Ej [C(tj + i)] 

+ B(ti, t i+1 )I JVt .> 7 v( ti )Ei [(M(t i+1 ) - K) + ] (54) 

5 Numerical results 

In this section, we present results obtained with BSLP based on calibration to iTraxx 
Europe Series 6 tranche quotes from March 15th, 2007. We look at two specifications of 
the model: 

A: The 0% ^-volatility case with mean reversion speed a = oo and transient contagion 
7 = 0. This is equivalent to the local intensity model presented in Sect. |2j 



B: A log-normal specification for the Y t process (19) with volatility a = 70%, mean 



reversion speed a = 30% and transient contagion 7 = 0. 
5.1 Calibration 

Calibration of BSLP to the tranche quotes was done using the 2D calibration described 
in Sect. 3J5 We present the quality of the calibration in Fig. [T] where we indicate mid as 



well as bid and ask levels. We see that in all cases the calibrated values lie well within 
the bid/ask spreads. 

5.2 Pricing 
5.2.1 Tranchlets 

Next we investigate the behavior of the two calibrated BSLP models A and B. First we 
look at implied tranchlet pricing, i.e. prices of thin CDO tranches on the iTraxx portfolio. 
In particular we look at the 10 year 1%-wide tranches up to 12% and plot the results in 

Fig. m 
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The first observation is that the tranchlet curve is close to linear on a logarithmic 
scale, resulting in a smooth interpolation and extrapolation to the equity region. 

The second observation is that the 0% and 75% volatility cases produce nearly iden- 
tical curves. This is a reassuring result that suggests that the local intensity model is 
indeed sufficient for tranche price interpolation. Note that this is not a completely trivial 
observation. Indeed, similarly to regular tranches, tranchlet prices are determined by a 
set of marginal loss distributions at a set of horizons. However, market tranche quotes do 
not uniquely define those marginal probabilities; they only provide integral constraints on 
them (i.e. the market is incomplete - see also the remarks at the end of section 3.6). When 
we calibrate BSLP with volatilities 0% and 75%, strictly speaking we end up with two 
different sets of marginal loss distributions, whose difference could in principle amplify 
when projected onto thin tranchlets. Thus our findings imply that volatility assumptions 
have a negligible impact on marginal loss distributions even for incomplete markets. 



5.2.2 Tranche Options 

The main aim of the BSLP model is to price exotic credit derivatives. Here we look at the 
tranche option as an example. We focus on 5y into 5y options, i.e. exercise date 20-Dec- 
2011 and maturity 20-Dec-2016. The tranches we consider here are the index (0% — 100%) 
and the equity (0% — 3%) tranche. 

Figs. [3] and [4] show the implied (Black) volatilities for these options for both model 
choices A and B. 

We start by discussing the index option (Fig. (13])). It can be seen clearly that for both 
models A and B there is a strong skew in the implied volatilities. 

It should be noted that even for model A which has 0% volatility for the portfolio 
intensity the implied volatilities are high (the ATM vol is around 22%). As we will see in 
the next section, this is because the model gives rise to a non-zero dispersion of forward 
spreads even in the absence of Y-volatility. 

One can also see that adding non-zero y-volatility results in a rotation of the volatility 
smile curve, i.e. volatilities increase for low strikes and decrease for high strikes. Again, 
this can be understood by looking at the behavior of the conditional forward spreads as 
discussed in section 15.2.31 

For the equity option (Fig. Q) we note that increasing Y- volatility results in an 
upward parallel shift of the smile curve. We can also see that the curve for the local 
intensity model (A) is not very smooth. This is a result of the fact that the local intensity 
model is fundamentally discrete, i.e. the dynamic default variable can only take on integer 
values. This means in turn that the distribution of conditional forwards is also discrete. 
As a result option prices are linear for strikes that lie between the conditional forward 
levels. 



5.2.3 Conditional Forwards 

A key quantity for assessing dynamics of a multi-period model is the forward spread 
distribution. Here we focus on the 5 year spreads, conditional on the default level at 
20-Dec-2011. 
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These spreads tell us how the model links spread and default dynamics. This can be 
compared to intuition on how risk will be priced in different (default) states of the world. 
Of course these forwards will also play a direct role in exotic pricing, e.g. a tranche option 
can be viewed as essentially an option on the forward spread. 

The forwards levels are plotted in Figs. [5] and [6] for 0% — 100% and a 3% — 6% tranche 
respectively. 

First we should note that 0% volatility (model A) does not imply variance for 
the distribution of conditional forwards. This leads to non-zero implied tranche option 
volatilities as described in the previous section. The reason for this is that even though 
there is no spread diffusion, we still have spread jumps determined by the contagion 
factors. 

One can also see a rapid increase in spread levels as defaults increase, which indicates 
that the model implies a high level of contagion. An interesting feature is that the increase 
in forward spreads is less pronounced for model B than A, i.e. introducing F-volatility 
reduces the level of contagion. 

This can be understood since we are replacing spread contagion with spread diffusion 
while calibrating to the same portfolio loss distribution. Another way to see this is to 
note that adding spread diffusion will decorrelate losses and spreads. 

This is an important effect, which is likely to have significant impact on exotic pricing. 



5.3 Risk 

In this section we explore the first order (spread) risk in the BSLP model. Since BSLP is 
driven by a stochastic F-process, it is natural to look at deltas with respect to moves in 
todays value Yq. The meaning of this procedure can be clarified using Eq.(21): if we keep 
contagion factors F(Nk, 0) constant and shift Yq, this is equivalent to shifting the index 
swap intensity A s (0), and thus can be used to specify the index delta. 

We note that, the same Eq.(21) also shows that shifting Yq while keeping contagion 
factors F(Nk, 0) constant is equivalent to a common rescaling of all contagion factors. 



5.3.1 Tranche Deltas 

First we look at how tranche prices behave as we shift the initial F-value, Yq. In particular, 
we can look at the hedge ratio to the index. 

Let us denote the mark-to-market of a tranche with low strike k, high strike I and 
coupon c as a function of Yq by: Ji4[k,i](c,Y ). The par spread of this tranche is S[k,i], 
defined by: 

M m (S m ,Y o ) = (55) 

For a given tranche with strikes k and I, we now define the tranche delta An^i as 
follows: 

A = M M (jjM ,Y + e)-M [kfi (Sjkfi , Yq) = M [kfi (S [kA ,Y + e) 
M " M m {S m ,Y + e)-M m (S m ,Y ) M m (S l0>1] , Y + e) 1 j 

For index tranches these deltas can be compared to the quoted spread deltas. This is 
done in Fig. [7j We see that deltas are reasonably close to the quoted values. BSLP deltas 
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tend to be lower in the equity region and higher for super senior tranches. We note also 
that adding y-volatility (model B) tends to improve the delta match. 

5.3.2 Option Deltas 

We now look at deltas of tranche options. Again we consider the impact of a shift in 
initial F-level and compute the option delta as the change in mark-to-market of the 
option divided by the change in mark-to-market of the forward tranche. In other words 
we are computing the hedge ratio of the tranche option with respect to a position in the 
underlying forward. 

We can compare this to the hedge ratio obtained by using the Black formula. In this 
case we first calculate the implied Black volatility of the option. After shifting Y we get 
new forward levels and risky annuities for the underlying tranche. Using these shifted 
forwards and annuities but keeping the implied volatility constant we can compute a new 
option price using the Black formula. This gives rise to a Black delta or hedge ratio. 

The results are given for a 0% — 3% equity tranche option and the 0% — 100% index 
option in Figs. [9] and [8] respectively. As before the exercise date of the options is 20-Dec- 
2011 and the maturity is 20-Dec-2016. 

We can see that BSLP and Black deltas are of the same order. The match is better 
for the equity option than the index option, especially for high strikes. We note also 
that Black and BSLP deltas are closer to each other for model specification B, i.e. with 
non-zero F-volatility. 

As in the case of the implied Black volatility smile we can see that the curve of option 
deltas as a function of strike is not very smooth for the 0% volatility case. Again this is 
due to the fundamental discreteness of the local intensity model. 

6 Summary 

BSLP is a model of the portfolio loss and loss intensity. It can be used for the pricing and 
risk management of vanilla as well exotic credit derivatives which depend on the portfolio 
loss. 

Because the model is low-dimensional and Markovian, efficient lattice or tree imple- 
mentations are possible. In addition to such numerical methods, analytical approaches 
such as the adiabatic approximation can be used. A possible extension to the model 
considers multi-factor generalizations of the driving stochastic process. 

BSLP explicitly models credit contagion and provides insight on how this impacts 
the forward loss dynamics. This is a key for the pricing and risk management of exotic 
portfolio credit derivatives. 

The model achieves a near perfect calibration to any set of portfolio tranche quotes 
due to a semi-parametric representation of the contagion factors. These quotes can either 
be obtained from the market for liquid indices or from an underlying bespoke tranche 
pricing model. 

Efficient lattice or tree pricing algorithms for the model have been developed, including 
fast calibration algorithms based on a ID projection of the full 2D model. Tree pricing 
is most suitable for pricing exotic credit portfolio derivatives that incorporate optionality 
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and/or weak loss-path dependence such as tranche options, forward tranches, leveraged 
super-senior etc.. However, Monte Carlo implementations of the model (not discussed in 
this paper) are also possible in order to price more strongly path dependent products. 

The calibrated model retains flexibility in adjusting the dynamics by choosing a con- 
crete specification of the underlying stochastic driving process. This provides some free- 
dom in the pricing of exotic derivatives. The local intensity limit of BSLP is a model in its 
own right that can be used to price non-standard tranches by arbitrage-free interpolation. 
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Calibrated values 
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1 2% 
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1 2% 


22% 
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0.40 


NA 


0.41 


0.34 


20-Dec-09 


22% 


1 00% 


NA 


0.10 


NA 


0.18 


0.12 


20-Dec-09 


0% 


1 00% 


1 1 .00 


12.50 


1 4.00 


12.39 


12.34 


20-Dec-1 1 


0% 


3% 


1 1 .75% 


11.88% 


1 2.00% 


11.88% 


11.88% 


20-Dec-1 1 


3% 


6% 


53.75 


54.50 


55.25 


54.50 


54.51 


20-Dec-1 1 


6% 


9% 


1 4.00 


14.75 


1 5.50 


14.75 


14.76 


20-Dec-1 1 


9% 


1 2% 


5.75 


6.25 


6.75 


6.25 


6.26 


20-Dec-1 1 


1 2% 


22% 


2.1 3 


2.50 


2.88 


2.52 


2.52 


20-Dec-11 


22% 


100% 


0.80 


1.05 


1.30 


1.18 


1.20 


20-Dec-11 


0% 


100% 
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0% 


3% 
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3% 


6% 
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6% 


9% 
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9% 
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5.50 


6.00 


6.50 


6.01 
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3% 
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6% 


9% 
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9% 
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40.00 
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4.60 


4.85 


4.65 


4.69 
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0% 
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45.00 
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Figure 1: Calibration results for BSLP. Note that quotes are given as spreads in bps 
except for the 5,7, and lOy equity tranches which are quoted as upfront in percent (with 
a coupon of 500bps). Calibration is done to mid quotes. Indicative bid/ask quotes are 
included for reference to assess calibration quality. Note that 3y bid/ask spreads are not 
available. 
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Tranchlet prices 
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Figure 2: Tranchlet (1%-wide) prices for models A and B on logarithmic scale. Maturity 
is the 20th December 2016. 
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Implied volatilities for 0%-100% tranche 

60% -i— 
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0% -I , 1 1 1 

1 2 3 4 5 

Strike as multiple of forward 

Figure 3: Implied volatilities for index (0% — 100% tranche) option. Exercise date is 
20-Dec-2011 and maturity is 20-Dec-2016. Strikes are given as a multiple of the forward 
value, which is at 69bps. 
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Implied volatilities for 0%-3% tranche 
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Figure 4: Implied volatilities for equity (0% — 3% tranche) option. Exercise date is 20- 
Dec-2011 and maturity is 20-Dec-2016. Strikes are given as a multiple of the forward 
value, which is at 2725bps. 
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Figure 5: Forward spreads on a logarithmic scale conditional on the default level for the 
index tranche for models A and B. Start date is 20-Dec-2011 and maturity is 20-Dec-2016. 
The unconditional forward level is ca. 69bps. 
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Fwd Spreads (bps) 
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Forward spreads conditional on defaults for 3%-6% tranche 
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Figure 6: Forward spreads on a logarithmic scale conditional on the default level for 
3% — 6% tranche for models A and B. Start date is 20-Dec-2011 and maturity is 20-Dec- 
2016. The unconditional forward level is 750bps. Note that the tranche is wiped out after 
11 defaults. 
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;ure 7: Tranche deltas computed with models A and B compared to iTraxx quoted 
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Option deltas for 0%-100% tranche 




Figure 8: Option deltas for the index tranche computed with models A and B compared 
to Black option deltas. Strikes are given as multiples of the forward level which is ca. 
69bps. 
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Option deltas for 0%-3% tranche 
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Figure 9: Option deltas for the 0% — 3% equity tranche computed with models A and B 
compared to Black option deltas. Strikes are given as multiples of the forward level which 
is ca. 2725bps. 



33 



Appendix A: Swap intensity and NtD intensity 

Here we formally derive the relation (pT| between the NtD intensities specifying the 



transition probabilities (17) in the full 2D model, and the stochastic index swap intensity 

The swap intensity is defined as the intensity that prices the index correctly as a single 
name CDS. More precisely, we want the following. Given a stochastic intensity X s (t) we 
denote the distribution, Q(t, T), of the first jump r of the associated counting process by: 

Q{t,T) = P[r>T\T t ] (A.l) 

This is all we need to price a CDS. To ensure that the value of this CDS corresponds to 
the value of the index calculated using the portfolio default process N t we need: 

Q[t , T) ,m^m (A . 2) 



The intensity corresponding to Q(t,T) is given by: 

dQ 



(t,T) (A.3) 



Using equation (A. 2) this gives: 



K[t)dt = m±*W - jVt (A .4) 



We can evaluate the conditional expectation in (A. 4) as follows: 

N-N t 



E[N t+dt \F t ] = ( N t + k)P[N t+dt = N t + k\N t ,Y t ,t] (A.5) 



k=0 



where P[N t +dt = N t + k\N t ,Y t ,t] is the probability of the transition Nt — > N t +dt = N t + k. 
As we only allow for at most one step transition in the infinitesimal time dt, the sum in 



(A.5) reduces to just two terms: 

E[N t+dt \F t ] = (N t + l)X NW {N t , Y t , t)dt + N t (l- \ NW {N t , Y t , t)dt) (A.6) 



Substituting this relation into ( A.4[ ), we arrive at the sought-after relation between the 
two intensities: 

A,W,n,t)= Af "° W ^"" (A.7) 

Note that X s (t) is also the average of the portfolio single name intensities Aj. This 
follows since in the absence of simultaneous defaults we know that the portfolio default 
intensity \NtD is just given by the sum of single name intensities. In other words: 

N-N t 



\m D = Yl ^ = (N-N t )X s (A.8) 
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Appendix B: BSLP in continuous time 

Here we present a continuous-time formulation of the BSLP model for a jump-diffusion 
specification of the driving V-process given in Eq.([T8|). (A continuous-time formulation 
with a discretized Y- variable is further analysed in Appendix C). In this appendix, we 
will use a more conventional notation (x, y) instead of (N, Y) for the dynamic variables 
of the BSLP model, with x denoting the defaulted fraction N t /N instead of the absolute 
default level N t . 

The general form of the infinitesimal jump-diffusion generator £ acting on the pdf 
P(z^ ,t\z, s) for n-dimensional random variables z 1 ', z G R n is as follows 

* • " ' d 2 P(z',t\z,s) 



£P(?,t\z,s) = s ) dP{ ^ S) +J2 g ) 



dzi ^— ' 2 y dzidzj 

+ W(z"\z,s) [P(Z,t\?',s) - P(z',t\z,s)} (B.l) 

while the adjoint operator £* readi 

n 8 - 1 <9 2 

i=l % i,j=l 1 1 

+ [W(?\z",t)P(?',t\z,s) -W(z"\z',t)P(z',t\z,s)] (B.2) 

z"^z 

Here W^lz, t) is the jump measure determining the jump probability in time dt together 
with the jump size distribution. Using standard conventions, we write 

W(2\z,t) = \(B,t)w(?-z\z,t) (B.3) 

where X(z, t) stands for the jump rate (intensity) and w(6z\z,t) is a pdf of the jump size 
5z given the initial point z. The forward and backward equations are 

We note that the generator can be written in the following form 

A = C + d (B.5) 

where £ an d C-i correspond to diffusion and jump parts of the generator, given by the 
first two terms and the last term in (B.l), respectively. 

23 Recall the definition of adjoint operator £*: 

JdxP l £[P 2 } = j dxP 2 C*[P 1 ] 
where P±(x), P 2 {x) are two arbitrary probability densities. 
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Let us now define the operators Co and C\ in our specific 2D setting with z = (x,y). 
In the BSLP model, the generator Co acts only on y. We therefore set 

d Id 2 
C P(x,y,t\x ,y ,s) = fx(y , s)—P(x,y,t\x ,y , s) + -cr 2 (y , s)—,P(x, y, t\x , yo, s) 



dy ' 



9yo 



(B.6) 

where the functions fi(y,t) and <r(y,t) are defined in the SDE (18). For the adjoint 
operators we obtain 



C* P(x,y,t\x ,y ,s) 



d Id 2 

[fi( y y,t)P( y x,y,t\x ,y ,s)} + -— [cr 2 ( y y,t)P(x,y,t\x ,y , s 



dy 



(B.7) 

The second generator C\ corresponds to the case where jumps x — > x + Ax are enabled, 
and in addition are accompanied by a related jump in y (which arises due to the dN t -term 
in (18)). The rate of this process is proportional to the product of the ^-variable and the 
fraction of surviving names, times the contagion factor F(x) = q(x)f(x), similarly to the 
treatment above: 

X{x,y,t) = y(l - x)q(x)f(x) 
and the jump size distribution is a product of two delta-functions: 



wi(x -x ,y- y \x , y , s) = 6(x-x - Ax)S(y - y - Ay) 
Using (B.9) and ( |B.8 ), we calculate the generator C\ 



(B.8) 
(B.9) 



dP(x,y,t\x ,y ,s) 



X(x Q , y , s) [P(x, y,t\x + Ax, y + Ay, s) 
P(x,y,t\x ,y ,s)} 



(B.10) 



and for the adjoint operator we have 



C[P{x,y,t\x ,y ,s) 



X(x - Ax, y - Ay, t)P(x - Ax, y - Ay, t\x , y , s) 
X(x,y,t)P(x,y,t\x ,y ,s) (B.ll) 



To summarize, Co in (B.5) collects terms in the generator C where Y t serves as the only 
dynamic variable while the default counting variable N t enters as a parameter. The second 
term C\ describes the part of the generator that contains both Y t and N t as dynamic 
variables. This special structure of the generator can be used to set up an adiabatic 
perturbation theory for the model (see Appendix D). 



Appendix C: Discretized BSLP in continuous time 

Here we discuss how the continuous-time formalism presented in Appendix B changes if we 
keep time continuous but assume that the Y" r variable can only take discrete values from 
some finite set. Such analysis may be of interest as in this case the model becomes that 
of a 2D Markov chain, allowing powerful methods of Markov chain theory to be used to 
build fast and accurate numerical approximations to the dynamics of the original model, 
get more insight into the model behavior, and build connections to previous research. 
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In particular, we will use this reformulation to compare our approach with the Davis-Lo 
model [12], as well as establish a link with a general and rich class of quasi-birth-death 
(QBD) processes. 

Unlike the setting of Appendix B, which uses a jump-diffusion framework, here we 
assume that the Y t -variable is discretized. For the general case of a two-dimensional 
process in the iVY-plane, the states are parameterized by two indices (i, n) (i for number 
of defaults, and n for Instate). Correspondingly, the 2D generator carries four indices 
instead of two. We'll use the notation A in \j m (t) for the matrix elements of the generator. 
In matrix notation, the generator can be viewed as a block matrix 



(C.l) 



where all matrices F^ % > have dimension M x M, with M being the dimension of 
the Y-space. The interpretation of these matrices is as follows. The matrix gives 
intensities of transitions between Y-states when there is no change of the L-state during 
the infinitesimal time step At, while the matrix provides intensities of joint events of a 
jump in the loss variable during the interval [t, t+dt] accompanied by a transition between 
Y-states. As usual, all off-diagonal elements should be positive, diagonal elements should 
be negative, and each row in A should sum up to zercp^j 

Given the 2D generator matrix A, the forward equation takes a block-matrix form 
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(C.2) 



where P(t,T) has matrix elements Pi n u m {t,T). 

The Davis-Lo model [12] is a particular case of a QBD credit process where the dy- 
namics of Y is locally deterministic. Matrix elements of A are given in this model by 2 x 2 
matrices 
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(C.3) 



Note that the form of F^ implies that whenever there is a default in a given time step dt, 
the probability for the hidden variable Y to stay in state "0" is zero, meaning it jumps 
to state "1" ("risky state") with (conditional) probability l 25 , or stays in state "1" with 



24 Such a block-matrix structure of the generator is characteristic of the so-called Quasi-Birth-Death 
(QBD) processes. In the terminology of QBD, the loss variable is the "level" while the Y-variable is 
the "phase" . Symbols "L" and "F" stand for "local" (without change of level) and "forward" (level is 
changed by one unit), respectively. 

25 Note that the clement — X(N — i) specifies the probability of the joint transition P[(i, 0) — > 
(i + 1, 1)] = A(7V — i)dt. On the other hand, the marginal transition probability P[i — ► i + 1] is also equal 
to A (TV — i). Therefore, the conditional probability P[0 — > l|z — ► i+ 1] = 1 as expected. 
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(conditional) probability 1 if it was already there at the beginning of the time step. The 
matrix is interpreted similarly. 

In the BSLP model, the dynamics of Y is Markov as opposed to being locally deter- 
ministic (as in the Davis-Lo model), thus matrices L^ 1 ' and will have very few, if any, 
zero elements, as in general all transition probabilities between different F-states will be 



non- vanishing. Note that the generator (C.I) can be identically re-written as follows: 
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where L® = L® + diag (F@l) , = diag (F®l) . Both matrices A and A x can 
be separately interpreted as generators, as in both of them off-diagonal elements are 
positive, diagonal ones are negative, and the row-wise sums are zeros. Note that they 
have a straightforward interpretation: generator A corresponds to the case when no 
L-transitions occur during the infinitesimal time step dt, while A± allows jumps in the L- 
process during the interval [t, t + dt], that are accompanied by possible transitions between 



y-states. Note that this decomposition is a discrete counterpart of (B.5). 



Appendix D: Adiabatic approximation in the contin- 
uous time BSLP model 

To get a tractable approximation to pricing index (or tranche) options, we assume that 
the characteristic time scales of changes in Y and iV {r spread and r loss , respectively) are 
significantly different. We expect this approximation to yield accurate results because 
spreads change daily while loss states change very infrequently. To acknowledge the 
presence of a hidden small parameter ~ T sprea d/Ti oss in the problem, we re- write the 



infinitesimal generator (B.5) of the BSLP model as follows: 



A=-C Q + d (D.l) 

e 

where the parameter e > is assumed to be small. Matrices Co/e and C\ correspond to 
parts of the generator responsible for the y-dependent part (where the default counting 
variable N t may enter as a parameter) A and a second term A\ that contains both iV t and 



Y t as dynamic variables. In particular, for the OU specification (20) with mean reversion 
parameter a—l/e and variance v 2 = a 2 /(2a), we obtain 

c = (0-yo)J^ + v 2 -f- 2 (D.2) 
dy dy 2 

Given these assumptions, we consider the adiabatic approximation (see e.g. [16] for 
financial applications of this method) to find the forward and backward dynamics in the 
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form of asymptotic series in powers of e. We restrict ourselves to a calculation of the 
leading order adiabatic approximation for the backward equation. 

Consider the backward equation for some function B(x,t) of initial values (x,t) (this 
function can be a conditional expectation or a transition probability viewed as a function 
of the initial variables): 

° D -CB = - (-Co + Ci) B (D.3) 



dt ' \e 

We look for a solution in the form 



B = B + eB x + e 2 B 2 + ... (D.4) 



Plugging this into (D.3) and equating like powers of e, we obtain the chain of equations 

CoB = 

C B l = — (D.5) 
dt 

C B 2 = — — — C\Bx 
dt 

As the generator C depends only on y, the solution to the first equation is 

B = B (x,t) (D.6) 

i.e. a function (unknown at this stage) of x and t only. To find this function, we multiply 
both sides of the second equation in (D.5) by the stationary state po(y,x) of the adjoint 
generator C , i.e. CoPo(y,x) = 0. The left hand side of this equation vanishes, hence the 
right hand side should vanish as well. This implies that B Q (x, t) should solve the equation 

0=(^ + £ 1 B > = ^? + (A>Bo (D.7) 

where 

(£1) = / dyp (y,x)Ci(y,x) (D.8) 



That is, Bq(x, t) satisfies the backward equation in the L-space where the fast F-dynamics 
is integrated out in the generator C\ with weight equal to the invariant distribution po(y, x) 
of Cq. Correction terms B\, B 2 , . . . can be found from further equations in (D.5). 

Let us consider the leading order adiabatic approximation formulas. Using (D.7) and 
(B.10), we obtain the following equation for the lowest order term B (x,t): 

= -X(x) (B (x + Ax, t) - B (x, t)) (D.9) 

at 

where 

\{x) = (1 - x)q(x)f(x)(y) = (1 - x)q(x)f(x) / dyyp (y,x) (D.10) 



Thus to the leading order in the adiabatic approximation, we find that the backward 
dynamics can be described as the ID pure default dynamics with a rescaled intensity. 
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If we further impose the condition that the intensity of this effective ID model should 
be equal to the intensity (1 — x)Xf(x) of the original default-only model, we obtain the 
following expression for the drift adjustment factor q(x): 



q(x) 
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(y)x Jyp (y,x)dy 



(D.11) 



One notes a similarity between(D.ll ) and (41). The difference between them is that 
while (41) contains the full time-dependent distrubution 7r(Y t ,t\N t ), it is the steady state 
distribution p (y,x) = ]xm t -^oo 7r(Y t , t\N t ) that appears in ( D.ll[ ). 

Corrections to the zero-order result Bo = Bo(x,t) can now be found from equations 



(D.5). In particular, after the right hand side of the second equation in (D.5) is "centered" 
(i.e. its expectation is zeroed by the choice of B ), we can write this equation as follows: 



C G B l = {{C l )-C l )B 



(D.12) 



This is the so-called Poisson equation that can be explicitly solved to find the first cor- 
rection Bi, as described in 
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